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ABSTRACT 

The solar system’s Oort cloud can be perturbed by the Galactic tide and by individual passing 
stars. These perturbations can inject Oort cloud objects into the inner parts of the solar sys¬ 
tem, where they may be observed as the long-period comets (periods longer than 200 years). 
Using dynamical simulations of the Oort cloud under the perturbing effects of the tide and 
61 known stellar encounters, we investigate the link between long-period comets and encoun¬ 
ters. We find that past encounters were responsible for injecting at least 5% of the currently 
known long-period comets. This is a lower limit due to the incompleteness of known encoun¬ 
ters. Although the Galactic tide seems to play the dominant role in producing the observed 
long-period comets, the non-uniform longitude distribution of the cometary perihelia sug¬ 
gests the existence of strong - but as yet unidentified - stellar encounters or other impulses. 
The strongest individual future and past encounters are probably HIP 89825 (Gliese 710) and 
HIP 14473, which contribute at most 8% and 6% to the total flux of long-period comets, re¬ 
spectively. Our results show that the strength of an encounter can be approximated well by a 
simple proxy, which will be convenient for quickly identifying significant encounters in large 
data sets. Our analysis also indicates a smaller population of the Oort cloud than is usually 
assumed, which would bring the mass of the solar nebula into line with planet formation 
theories. 

Key words: stars — Oort Cloud — celestial mechanics — solar neighbourhood — Galaxy: 
kinematics and dynamics 


1 INTRODUCTION 

Planetary craters reveal an extensive history of bombardments 
caused by asteroids and comets. Large impacts on the Earth 
may have trigge red global catastrophic events such as the K-T 
mass extinction ^Alvarez et al J1 1 98(1 : ISchulte et alJl201(ih . To un¬ 
derstand such events and to predict the risk of future impacts 
it is necessary to study the past and current impact rate in the 
solar system . One approach to this is to date impact craters 
dlvano \][20Q2j;[Bottke e t al]2005l ). although many terrestrial craters 
have large age uncertainti es an d very few ages exist for extrater¬ 
restri al craters dwi lhelms. McCauley & Traskl 1 19871 : ICuller et ail 
l200d : Ijourdan. Reimold & Deu tschll2012h . Another approach is to 
dynamically simulate the motions of asteroids and comets and to 
calculate their impact rates. For example, the late heavy bombard¬ 
ment c an be explained by t he migration of gian t planets in the Nice 
model d Gomes et alJl2005l : lTsiganis et alj|2005h . 


cloud (OC), which was proposed by lOor 31195(3) to explain the 
large semi-major axes (> 20, 000 au) of long-period comets (those 
with periods larger than 200 yr). The angular momenta of comets 
in the OC are changed by the gravitational perturbations from 
the Galactic tide and stellar encounters. When the angular mo¬ 
mentum or perihelion of an OC comet becomes small enough, it 
would be observed and classified as a long-period comet (LPC). 
By integrating the orbit of an LPC back in time we may be able 
to identify the stellar encounter responsible, i Dvbczvhskil 1 20021 : 
iRickman et alJl2012l : iFeng & Bailer-Jone^ 12014) . In general, dy¬ 
namical simulations of the formation and evolution of the OC are 
crucial for improving our understanding of the origin of Sedna-like 


bodies (Brown, Trujillo & Rabinowitzll2004l: Morbidelli & Levisonl 


200j^jKei^onjfcRromley| 20041 : iKaib. Roskar & Ouinnl 201 ll : 


Trujillo & Sheppard 2014), the birth environment of the Sun 
dLevison et alJl201 0: Brasser et al. 2012), and the cometary impact 
rate on the Earth dFeng & Bailer-Jone5l2014lf . 


In addition to the analysis of minor bodies in the inner so¬ 
lar system, many previous studies have examined the evolution 
of the outer solar system and the planetary impact rate asso ciated 
with it teverharil967l : IWeissmanl20'07l : lKaib & Quinn|[2009h . This 
may be modulated by the solar orbit in the Ga laxy, as this changes 
the stellar density in t he sol a r neighbourhood icia rdner et ali i 20 111 : 
IFeng & Bailer-Jonesl 120131 : IFeng & Bailer-Jones] 2014ll . Specifi¬ 
cally, the outer part of the solar system is occupied by the Oort 


However, few studies have systematically investigated the in¬ 
fluence of observed stellar encounters on the OC and thus on 
the flux of LPCs. Most studies have only estimated the influence 
of a few close encounters or have used empirical models, rather 
than through dynamical simulations. In this work we try to bet¬ 
ter quantify the influence of encounters by using dynamical sim¬ 
ulations of the OC to model the cumulative effect of many close 
encounters. We collect known close and strong encounters from 
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iBailer-Jonesl J20 1 5t) : IDvbczvnski & Berskil j20l5t) : lMamaiek et alj 
d2015l) . We then simulate the OC comets by taking into account the 
gravitational perturbations from these encounters and the Galac¬ 
tic tide. To investigate the role of each perturber, we conduct sim¬ 
ulations of the OC under perturbations from the tide only, from 
encounters o nly, and from both (which we call “combined simu¬ 
lations”; see iFeng & Bailer-Jonesll2014l for details). We use sub¬ 
scripts t, e, c to refer to these types of simulations respectively. 
Using these results we measure the effect of each encounter us¬ 
ing various metrics. By comparing the results of the three types 
of simulations we investigate the nature of non-linear interactions 
between the tide and encounters, the so-called“synerg y effect” 
t lHeisler, Tremaine & Alcockll 19871 : [Rickman et ~al]I2008B . Finally, 
we compare the angular distributions of the perihelia of LPCs and 
encounters to assess the influence of the solar motion on the current 
LPC flux. 


2 DATA 


We collected the data of stellar en counters from three sources : 
iBailer-Jo nes! d2015l) thereafter BJ15). lDvbczvnski & Berskil J20 1 5l) 
and Mamaiek et alj J2015h . Following Bailer-Jonesl 1 2015h . we 
use the term “object” to refer to each encountering star in our 
catalogue. A specific star may appear more than once but with 
different data, thus leading to a different object. As described 
in BJ15, we use GCS (G), RAVE (R), Pulkovo (P) to denote 
the catalog of objects with radial velocity obtained by cross¬ 
matching the revised Hipparcos catalog tvan^eeuwen 2007) with 
the Geneva-Copenhagen surv ey data bv dC asagrand e et alj|20llh . 
with the Pulkovo catalogue dGontcharovI 20061) . and with Rave- 
DR4 dKordonatis et alj|2 013 ), respectively. We als o use data from 
the XHIP (X) catalo g dA nderson & Francis] 2012h and the XF1IP 
catalog adapted by ( IDvbczvnski & Berskil 20151) (D). We also 


include WISE J072003 20-084651.2 (W0720) wi th data from 
dScholz. Bihain & Stormll20l4l ; lMamaiek et al1l2015l) (M). In sum¬ 
mary, the catalogs of G, R, P and X are from BJ15, and catalog D 
and M are from their corresponding references. 

We first select all objects with mean perihelia d p h < 2 pc from 
the catalog provided by BJ15. This gives 65 objects. However, more 
distant encounters can also have significant perturbing effect. We 
therefore select additional encounters using a proxy, g, for the size 
of the perturbation, defined as 


= M en c 

9 d lhVph 


7 

d ph 


( 1 ) 


where M enc is the mass of the encountering star in solar masses, 
d p h is its perihelion distance in pc, and u D h is its velocity at per¬ 
ihelion in the heliocentric rest frame (Rickmaij tl976 ) in km/s. 
7 (which is defined by the above equation) is linearly propor¬ 
tional to the impulse of the Sun gained from the encounter, and 
was used as a proxy for en counter strength in our previous work 
dFeng & Bailer-Jonesll20l4) (hereafter FBJ14). We will see in sec¬ 
tion |4j] that g is actually a better proxy, which is why we use it 
here. 

The values of g for the 65 selected objects ranges from 10 -4 to 
3. Of these, 34 objects have g < O.Ol(M 0 pc _ 2 km _ 1 s), and these 
can only inject about one in a million LPCs from the OC comets 
(see section ldTI and Table!]}. To select encounters with at least this 
perturbation strength but which have perihelia beyond 2 pc, we ad¬ 
ditionally select all encounters with g > O.Ol(M 0 pc 2 km 1 s), 


yielding another 25 objects from BJ15’s catalog. We then also se¬ 
lect the 1 2 new encountering ob j ects w ith perihelia less than 2 pc 
found by iDvbczvn ski & B erskil l2015i), as well as the close en¬ 
counter W0720 from Mamaiek et alj d2015h . The masses of these 
objects are either from the literature or estimated according to their 
stellar types based on Table 1 of lRickman et al.1 d2008h . Since many 
sources may be unresolved binaries, the masses may have uncer¬ 
tainties as large as 200%, especially for the estimated masses. We 
are left with a sample of 103 objects, of which 90 are from BJ15 
and 13 are from other catalogs. We will do full integrations for all 
selected sources. 

For each object, we resample the observed data according to 
the observational errors (as described in BJ15) to create 100 dif¬ 
ferent initial conditions. We integrate the orbits of all samples, and 
select the trajectory which gives a perihelion closest to the mean 
perihelion given in the catalogs Q. This trajectory is more repre¬ 
sentative of the mean perihelion than using the nominal observed 
data, due to the nonlinear propagation of errors (see BJ15). The 
time (f p h), velocity (i> p h), and distance (d p a) of each object at its 
perihelion are calculated for the selected trajectory. The resulting 
data are shown in Table|T| 

In what follows we study the individual and collective influ¬ 
ence of the 103 encountering objects on the OC. For this collective 
assessment we select a unique object for each star to generate a 
sequence of encounters such that each star occurs only once. We 
select objects from the catalogs in the following decreasing order 
of priority (set by our assessment of their reliability): G, P, X, M, D, 
and R. We call this set of object the basic sequence. Some objects 
in this sequence have dubious radial velocities, proper motions, or 
parallaxes. B y excluding these objects acc ording to the comments 
in BJ15 and IDvbczvnski & Berskil J2015I) . we define a subset of 
the basic sequence which we call the reliable sequence. Encoun¬ 
ters are labelled “b” or “r” in Table [Tjaccording to which sequence 
they belong to (all “r” are also in the basic sequence), while those 
labelled“n” are not in either sequence. 


3 METHOD 

3.1 Simulation scheme 


Following the simulation scheme in FBJ14, we conduct three sets 
of simulations to study the influence of stellar encounters on the 
OC: tide-only, encounter-only, and c ombined encounter-tide sim¬ 
ulations using the AM USE s oftware i Portegies Zwart et al.ll200 c } : 
IPelupessv. Janes & Porteeies Zwartl 120121) . We use the Galaxy 
model and initial conditions of the Sun as in BJ15. The initial 
conditions of the OC comets are from FBJ 14, which were gen¬ 
erated from a model adapted from iDones et alj (i 2004 a)’s results, 
and which we call the DLDW model. In this we used an inner 
boundary for the smallest perihelia of comets in order to avoid 
the influence of planets on the OC. However, this results in a non¬ 
equilibrium initial distribution. To overcome this problem, we use 
tide and encounter-like perturbations (see FBJ 14) over 100 Myr to 
move the OC closer to the kind of steady state it would have around 
the present time. We then take the coordinates of 10 million comets 
which remain in the OC at the end of this as the initial conditions 
for the following simulations. We will integrate the orbit of each 


1 The mean and median of the distributions of perihelion time, distance and 
velocity are very similar: see BJ15 and its supplementary information on 
http://www.mpia.de/homes/calj/stellar_encounters/stellar_encount 
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Table 1. Data on the 45 most perturbing encounters according to the number of encounter-induced LPCs, N e (and sorted by this value). This value was 
calculated from the encounter-only simulations using one million OC comets. The objects are denoted using Hipparcos IDs except for W0720, which is 
denoted by “0720”. The penultimate column ("seq”) denotes the sequence which an object belongs to: “r” denotes objects in the reliable sequence, a subset 
of the basic sequence with reliable data; “b” means objects in the basic sequence but not in the reliable sequence; and “n” indicates objects in no sequence. 
The reference to the stellar masses is given in the last column: “(l)-estimated” means the mass was estimated according to the stellar type. The full list of 103 
encounter objects is available online. 


HIPID 

cat 

A^enc 

M 0 

iph 

Myr 

^ph 

km s —1 

d p h g x 100 

pc M0pc 2 km 1 

N e 

S 

/max 

XlO 3 

AT max 

Myr 

fo 

xlO 4 

seq 

ref 

85605 

X 

0.7 

0.34 

20.98 

0.103 

312.0 

539 

555.3 

0.80 

0.0 

b 

(1) 

89825 

X 

0.6 

1.42 

13.99 

0.267 

60.4 

187 

154.6 

1.19 

0.0 

n 

( 2 ) 

63721 

p 

1.7 

0.14 

36.09 

0.273 

63.2 

159 

137.5 

1.30 

0.0 

b 

(1) 

89825 

p 

0.6 

1.30 

14.76 

0.363 

30.8 

115 

82.0 

1.56 

0.0 

r 

( 2 ) 

14473 

D 

1.1 

-3.52 

36.16 

0.282 

38.4 

67 

58.4 

1.79 

362.1 

r 

(1) 

14576 

X 

6.0 

-8.48 

3.28 

4.059 

11.2 

38 

21.5 

1.43 

38.0 

n 

( 3 ) 

103738 

P 

2.5 

-3.71 

18.09 

0.825 

20.3 

31 

14.9 

9.33 

61.0 

r 

(4) 

14576 

P 

6.0 

-4.30 

6.48 

3.145 

9.4 

31 

18.0 

1.34 

114.4 

r 

( 3 ) 

85661 

P 

1.7 

1.92 

46.39 

0.594 

10.4 

22 

9.7 

5.10 

0.0 

r 

(1) 

84263 

D 

1.0 

-6.29 

11.92 

1.201 

6.0 

14 

7.6 

0.80 

36.5 

r 

(1) 

27288 

X 

3.2 

-0.86 

24.76 

1.300 

7.6 

12 

7.8 

5.14 

45.4 

n 

(1) 

27288 

P 

3.2 

-0.86 

24.50 

1.300 

7.7 

12 

7.8 

5.07 

45.4 

r 

(1) 

26744 

D 

1.3 

12.64 

4.83 

1.714 

9.2 

9 

2.2 

9.31 

0.0 

n 

(1) 

25001 

P 

3.1 

4.92 

17.00 

2.185 

3.8 

8 

6.4 

13.86 

0.0 

r 

(1) 

71683 

G 

1.1 

0.03 

34.14 

0.902 

4.0 

7 

5.1 

0.84 

0.0 

r 

(5) 

71683 

P 

1.1 

0.03 

34.29 

0.910 

3.9 

7 

5.1 

0.84 

0.0 

n 

(5) 

71683 

X 

1.1 

0.03 

31.84 

1.001 

3.4 

7 

5.1 

0.84 

0.0 

n 

(5) 

87052 

D 

3.2 

5.75 

35.68 

1.812 

2.7 

7 

5.1 

0.88 

0.0 

r 

(1) 

71681 

R 

0.9 

0.03 

29.67 

0.912 

3.6 

6 

5.1 

0.84 

0.0 

n 

(5) 

25001 

D 

3.1 

3.97 

22.83 

2.069 

3.2 

6 

5.7 

14.47 

0.0 

n 

( 1 ) 

30344 

P 

0.8 

-1.94 

14.71 

1.394 

2.7 

5 

3.7 

17.66 

7.5 

n 

( 1 ) 

71681 

X 

0.9 

0.03 

30.32 

1.062 

2.6 

4 

3.7 

19.55 

0.0 

r 

(5) 

30344 

G 

0.8 

-1.98 

14.25 

1.415 

2.7 

4 

3.7 

17.67 

6.5 

r 

( 1 ) 

26624 

X 

2.1 

-1.84 

22.89 

1.939 

2.4 

4 

3.7 

17.78 

8.8 

n 

( 1 ) 

26624 

P 

2.1 

- 1.88 

22.06 

1.984 

2.4 

4 

3.7 

17.74 

7.9 

r 

( 1 ) 

25001 

X 

3.1 

3.64 

23.47 

2.226 

2.7 

4 

3.5 

3.45 

0.0 

n 

( 1 ) 

42525 

X 

1.3 

-0.26 

59.99 

0.955 

2.4 

3 

2.1 

6.26 

20.3 

r 

( 1 ) 

94512 

X 

3.2 

3.77 

30.21 

1.829 

3.2 

3 

2.1 

9.09 

0.0 

n 

( 1 ) 

93506 

X 

5.3 

- 1.01 

25.90 

3.323 

1.8 

3 

4.1 

0.45 

35.6 

n 

( 6 ) 

93506 

P 

5.3 

- 1.02 

26.03 

3.329 

1.8 

3 

4.1 

0.43 

34.7 

r 

( 6 ) 

23415 

D 

3.2 

4.22 

30.90 

1.650 

3.8 

3 

2.1 

0.16 

0.0 

r 

( 1 ) 

91012 

R 

1.3 

0.26 

349.23 

0.486 

1.6 

2 

4.2 

5.90 

0.0 

b 

( 1 ) 

30344 

X 

0.8 

-1.55 

18.86 

1.160 

3.1 

2 

2.1 

17.31 

0.0 

n 

( 1 ) 

38965 

X 

1.3 

-1.23 

54.83 

1.330 

1.3 

2 

2.1 

0.02 

10.0 

r 

( 1 ) 

25240 

P 

1.5 

-1.06 

54.60 

1.623 

1.1 

2 

2.1 

0.44 

17.2 

n 

( 1 ) 

25240 

G 

1.5 

-0.99 

53.82 

1.665 

1.0 

2 

2.1 

0.50 

18.3 

r 

( 1 ) 

94512 

P 

3.2 

3.51 

29.59 

2.241 

2.2 

2 

2.1 

14.40 

0.0 

r 

( 1 ) 

32349 

P 

3.0 

0.04 

17.54 

2.505 

2.7 

2 

2.1 

0.48 

0.0 

r 

(7) 

32349 

X 

3.0 

0.05 

17.55 

2.498 

2.7 

2 

2.1 

0.48 

0.0 

n 

(7) 

26744 

X 

1.3 

12.98 

5.19 

2.630 

3.6 

2 

1.1 

-2.90 

0.0 

n 

( 1 ) 

26744 

P 

1.3 

16.13 

3.94 

2.990 

3.7 

2 

1.1 

-6.05 

0.0 

n 

(1) 

26744 

G 

1.3 

13.16 

4.74 

2.969 

3.1 

2 

1.1 

-3.08 

0.0 

r 

( 1 ) 

18907 

X 

2.2 

5.00 

6.93 

5.107 

1.2 

2 

3.9 

1.14 

0.0 

n 

(8) 

18907 

P 

2.2 

4.96 

7.05 

5.698 

1.0 

2 

3.9 

1.18 

0.0 

r 

(8) 

W0720 

M 

0.1 

-0.07 

83.58 

0.252 

2.8 

1 

5.2 

0.44 

52.3 

r 

(9) 


(l)-esti mated, (2)j Garcfa-Sanchez et al.Hl99^1. (3)jRhee et al J J2007I). (4)jBailer-Jonesl J 2 1 ) 1 di (5Y| Thevenin et al.| |2 l)02fi . 

(6t lDe Rosa et al.l j2012h . (71- lLiebert et al.l i2005l) . (8) jMalagnini & Morossiil990h . (9) jMamaiek et al .11 20151) 


individual OC comet under the influence of the gravitational per¬ 
turbations. 

In the following we apply the tide-only, encounter-only, and 
combined simulations to the time span covering each individual 
stellar encounter as well as the two sequences of encounters. For the 
study of individual encounters we simulate one million OC parti¬ 
cles over a time span of 10 Myr or 20 Myr from now, either into the 
future or the past according to when the perihelion occurs. Specif¬ 


ically, if |f p h| > 10 Myr we adopt a 20 Myr time span, otherwise 
we use a 10 Myr. In the case of encounter sequences, we simulate 
10 million OC comets over a time span of 10 Myr. Encounters with 
perihelion times out of this window are naturally excluded. 
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tph tmax 



Figure 1. The distribution of the flux ratio, /, over time for HIP 14576 (Al¬ 
gol) based on backward simulations of 1 million OC comets over lOMyr. 
This distribution is obtained from a kernel density estimate over the injec¬ 
tion times of the individual LPCs (shown in grey lines at the bottom) using 
a Gaussian kernel with a scale length of 1 Myr. The green and red lines in¬ 
dicate t = 0 Myr and t = t p (encounter perihelion) respectively. The blue 
line marks the time when f reaches its maximum. 

3.2 Analysis of LPC flux 

Having integrated the comet orbits in the above simulations, we 
then calculate the resulting LPC flux. To do this we use two defini¬ 
tions of an LPC. 

First, an OC comet is counted as an LPC when its perihelion 
is less than 15 au. Within this so-called “loss cone” (LC), plane¬ 
tary perturbations can st rongly influence the orbit of the comet and 
either capture or eject it ( IWiegert & Tremaine! 19991) . Although ob- 
served LPCs have perihelia less than 5 au (the “observable zone”), 
we found in earlier work (FBJ14) that the flux of comets injected 
into the LC is approximately linearly proportional to that injected 
into the observable zone. The injection time of an LPC is when its 
pericentric distance is equal to 15 au, r = q\ c . This injection time 
is calculated by extrapolating the trajectory of the comet from the 
place where it gains the impulse to its perihelion, assuming a Kep- 
lerian motion. 

In addition to LPCs in the LC (or LC LPCs), there is another 
type of LPC, “dynamically new LPCs”. These are comets which 
have enough energy to penetrate the Jupiter-Saturn barrier and en- 
ter th e observabl e zone within one revolution dHeisler & Tremaine! 
1 1986h . Following lDvbczvnski & Krolikowskal J2015h . we count an 
LPC as a new LPC when its perihelion decreases from a value 
larger than 15 au to a value less than 5 au within one revolution. 
Since planetary perturbations probably cannot change a new LPC’s 
perihelion significantly over one revolution, our neglect of plane¬ 
tary perturbations in the simulations does not significantly diminish 
the accuracy of the computed LPC flux. 

With these definition of LPCs we calculate the injection time 
of LPCs generated in simulations of the OC. We use kernel den¬ 
sity estimation to convert these discrete times into a continuous 
LPC flux as a function of time, F e (f). From this we can estimate 
the maximum and the current LPC flux induced by individual en¬ 
counters. To compare this with the effect of the Galactic tide, we 
estimate the strength of an encounter using the ratio of encounter- 
induced to tide-induced LPC flux. 

Since the tidal force imposed on the OC from the Galaxy 


hardly changes over a few million years (FBJ15), the LPC flux 
caused by the tide, Ft(t), can be considered as constant, i.e. 
Ft(t) = Ft over our integration time span of 10 Myr. To quan¬ 
tify the influence of encounters on the Oort cloud, we define the 
flux ratio as 

/(f) = F e (t)/F t , (2) 

where Ft is the average flux over the integration time span. The 
encounter-induced flux, F e (t), is derived from the encounter times 
using a Gaussian kernel with a scale length of 1 Myr. Since the tide 
dominates the production of cur rent LPCs dMatese & Whitmirel 
2011; Dvbczvnski & Berski 2015), the flux ratio is also a good ap¬ 
proximation of the contribution of encounters to the total LPC flux. 

To characterize the perturbations on the OC from encounters, 
we further define the maximum flux ratio as 

/max = Te(fmax)/ F t , (3) 

where F e ( f ma .x) is the maximum flux of the encounter-induced 
LPCs, and f m ax is the corresponding time. If the encounter-induced 
LPC flux is small, /max may be sensitive to the scale length used to 
smooth the temporal distribution of LPCs. We should exercise cau¬ 
tion in interpreting /max when the number of encounter-injected 
LPCs, N e , is less than 10. From this we also define the time delay, 
AT m ax, between the time when encounter-induced flux reaches its 
maximum and the perihelion time, f p h. 

To estimate the influence of encounters on the current LPC 
flux, we further define the current flux ratio as 

fo = F e (t = 0)/ F t . (4) 

Since an encounter can only significantly influence the OC at its 
perihelion - very slow encounters being rare exceptions - future 
encounters generally do not influence the current LPC flux. Thus fo 
is forced to be zero for future encounters. The calculation of /max, 
ATmax and fo is illustrated by the example of Algol in Figure Q] 


4 RESULTS 

4.1 LPCs induced by individual encounters 

We conduct tide-only, encounter-only, and combined simulations 
over 10 Myr or 20 Myr from now for each encounter for one million 
OC comets generated from the DLDW model. We count the result¬ 
ing number of LPCs as a function of time, analyze the influence of 
encounters using the statistics introduced in section [3~2l and show 
the results in TableQ] Among the 103 objects in our encounter sam¬ 
ple, we only show the quantities of the 44 strongest objects which 
induce at least 2 LPCs in the simulations. We also show the results 
for W0720, which/s one of the closest stellar encounters found so 
far jMamaiek et all 20 lH) . 

In table [T| there are 12 objects which induce more than 10 
LPCs. Of these, seven induce LPCs which have perihelion times 
around t = 0 Myr. We select all the objects in the reliable sequence 
and add their current flux ratio, fo, and find that these objects in¬ 
duce about 5% of the current LPC flux. Although the first three 
objects have stron g influence on the Oort cloud , their data are un¬ 
reliable (BJ15 and lDvbczvnski & Berskj||2015r) . so we do not dis¬ 
cuss their results further. We find that the encounter which has reli¬ 
able data and induces the most LPCs is HIP 89825 (GL 710). This 
object is also confirm ed as the strongest near-future encounter b y 
many previous studies dGarcfa-Sanchez et alj !999l ; lBobvlevl201(]| ). 
However, it can induce an LPC flux which is only 8% of the tide- 
induced flux. In contrast, dMatese & Lissaueill2002l) has found GL 
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710 to induce 1.4 times as many LPCs as the tide. This discrepancy 
is understandable because we calculate the time when comets ac¬ 
tually reach their perihelia, t(r = q\ c ), while dMatese & Lissaueil 
12002 !) use the time when comets achieve the necessary impulse to 
(later) reach a small perihelion, t(q = qi c ). Another way of seeing 
this is that we calculate the flux by effectively averaging the num¬ 
ber of LPCs over the typical orbital period of an LPC (about 2 Myr) 
while they effectively average over the encounter time scale (about 
0.1 Myr). 

Since HIP 856 05 and 63721 do not have re liable astrometric 
solutions (BJ15 and Dybczynski & Berski 2015), we instead con¬ 
sider GL 710 and 85661 as the strongest two future encounters. The 
strongest past encounter is HIP 14473, contributing about 3.6% to 
the current LPC flux, and at most 6 % of the total LPC flux. We 
note, however, that this star has larger proper motion uncertainties 
fovbczvhski & Berskil201 if) . We find that the other two strong past 
encounters are HIP 103738 (gamma Microscopii) and HIP 14576 
(Algol), both of which are massive and relatively slow. The radial 
velocity of Algol is very close to zero with a standard error sim¬ 
ilar to its magnitude, leading to a large uncertainty in the perihe¬ 
lion time. Our results are largely consistent with the findings of 
iFouchard et all J201 ll) . namely that massive stars can have strong 
effects on the Oort cloud without having to make a very close en¬ 
counter. 

To illustrate further the role of encounters in perturbing the OC 
comets, we take HIP 89825 (GL 710) as an example. The upper 
panel of Figure [2] shows the temporal distribution of the induced 
LPCs and reveals a long-term influence of HIP 89825 on the OC. 
In the lower panel we see that HIP 89825 can change the perihelion 
of an LPC instantaneously and thereby allow the tide to gradually 
reduce the perihelion until the comet is captured by the planets. 

It is interesting also to examine the imprints of encounters on 
the distributions of the orbital elements of LPCs. Figure[3]compares 
the orbital elements of comets in the initial OC with those of the 
comets generated in tide-only, encounter-only, and combined simu¬ 
lations for HIP 103738 (Gamma Microscopii). (For the encounter- 
only simulations we increased the samples size to 10 million OC 
comets in order to improve the statistical significance.) The orbital 
elements of LPCs are calculated when they are injected into the LC; 
this is appropriate because most orbital elements do not change sig¬ 
nificantly over the course of a single orbit. 

We performed this analysis for other encounters too. Overall 
we find that encounters tend to push comets with large semi-major 
axes into the inner solar system (i.e. to become LPCs), because 
the outer part of the OC can experience stronger tidal forces from 
weak encounters (with large d p h). It is necessary that all LPCs have 
large eccentricity and small perihelia in order to be LPCs. How¬ 
ever, while the encounter-induced LPCs can have perihelia much 
less than 15 au, tide-induced LPCs have perihelia which concen¬ 
trate around the LC limit (15 au). This highlights the important role 
of encounters for injecting comets into the observable zone. 

In addition, tide and encounter perturbations have different 
influences on the other orbital elements. The tide tends to inject 
comets with high inclination while encounters do not. The LPCs 
induced by the tide generally show a W-shaped distribution in the 
longitude of the ascending node, H. Finally, the argument of per¬ 
ihelion of tide-induced LPCs concentrate around to = 50 deg and 
to = 240 deg while the distribution for LPCs induced by Gamma 
Microscopii is concentrated around to = 100 deg. These non- 
uniform distributions are caused by the anisotropic perturbations 
imposed on the OC from the tide and encounters, which we will 
further discuss in section l4~4l 
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Figure 2. Upper panel: the distribution of the injection time of LPCs in¬ 
duced by HIP 89825 (GL 710) for a simulation of 10 million OC particles. 
Lower panel: the variation of the perihelion of an example LPC under the 
perturbations from GL 710 (dashed black line), from the tide (thin black 
line) and from both together (thick black line). The red and green dotted 
lines show the upper limits of the loss cone and observable zone respec¬ 
tively. The blue vertical line shows the perihelion time of GL 710 in both 
panels. 


4.2 Comparison of perturbation strength proxies 

In section [2] we used the proxy, g, to select encounters. Here we 
compare the two proxies, g and 7, with the actual number of LPCs 
induced, N e . Figure [4] shows the log-log plot of the relation be¬ 
tween each proxy and N e . We fit linear functions of 7 and g with¬ 
out intercept to N e , i.e. N e = ay and N e = /3g, where a and 
/3 are unknown and to be fit for. The g proxy fits the number of 
encounter-induced LPCs with y 2 = 12.6 compared to y 2 = 47.1 
for 7, indicating that g is a better proxy. 

We also vary the index of d p h to find the best function form of 
N e which is 

N e = ?)Menc/(d p hPph) , ( 5 ) 

where g and S are unknown. The best fit of the data gives an ex¬ 
ponent of 1.92 (5 = 1.92) for d p h (Eqn.Q}. which also suggests 
that the g proxy is better than the 7 proxy. However, this analysis 
only includes observed stellar encounters, which naturally includes 
very few rare, strong encounters. It also omits very weak encoun¬ 
ters, which could not have perturbed any of the 1 million particles 
in our simulations. Thus 7 could still be a good prox y for v ery weak 
or very str ong stellar encounters, as is s uggested bv lFouchard et alJ 
J201 if) and lFeng & Bailer-Jonesl d2014l) . 
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Figure 3. The distributions of orbital elements of initial OC comets (black lines) and comets that are injected into the LC in three types of simulations for HIP 
103738 (other lines). The former are orbital elements of OC comets before any tidal perturbations while the latter are the ones calculated when comets are 
injected to become LPCs. The legend in the top right panel defines the lines for all panels: comets generated from the OC, and LPCs generated in combined, 
tide-only, and encounter-only simulations. In clockwise order from the top left, the panels show the: semi-major axis a,; eccentricity e; perihelion q; inclination 
I ; longitude of ascending node 9; argument of perihelion ui. Note that the distribution of cos I is uniform if the LPCs are isotropic. The number of comets in the 
OC and the LPCs generated in the combined, tide-only, and encounter-only simulations are 10®, 1892, 1912, and 444 respectively. To display the differences 
between encounter-induced and tide-induced LPCs better, we plot log densities over log |(] ( 1 — e) (to better distinguish between very high eccentricities) and 
l°glo(<j/au), rather than e and q. The OC model used in the encounter-only simulations contains 10 million comets. The range of orbital elements are divided 
into 100 and 20 bins for the original OC comets and LC LPCs, respectively. The strong fluctuations of the red lines around log 10 (l — e) = —3 and for 
l°glo(g/au) < 0.8 are caused by sample noise. 


4.3 LPCs induced by encounter sequences 

In addition to individual encounters, we combine encounters to 
form encounter sequences and simulate their influence on the OC 
over the past and future 10 Myr. To effectively account for the effect 
of weak encounters, we simulate with a sample of 10 million OC 
particles. We apply the tide-only, encounter-only, and combined 
simulations to integrate the motions of OC particles for both the 
reliable and the basic encounter sequences (see section^. The re¬ 
sulting number of LPCs is shown in Table [2] We see that over the 
past/future 10 Myr, the reliable sequence of encounters contributes 
about 8% to the total LPC flux (column N e /Nt). The basic se¬ 
quence contributes about 38% to the total future flux, due primarily 
to the strong perturbations from HIP 85605. Although HIP 85605 
has dubious astrometry (BJ15), we nonetheless include it in the 
basic sequence, not least to get an idea of the influence on the OC 
of possible undiscovered strong encounters. The results for the reli¬ 
able sequence should be used to see the effect of known encounters, 
although this sequence is almost certainly incomplete. 

Table [2] also shows that encounters seem less able to generate 
dynamically new comets by themselves (column n e /nt ). This is 


because they cannot change the perihelion sufficiently over a short 
interaction period (typically a few thousand years). The tide, in 
contrast, can gradually reduce the perihelia of comets by a large 
amount. Assuming there are 1 0 12 comets in the OC dWeissmanl 
1 199CL fl996: Pon es et al]|2004bl) . the averaged rate of dynamically 
new LPCs (with perihelia less than 5 au) is about 70 per year ac¬ 
cording to our simulations. For combined simulations for the reli¬ 
able sequence, there are about 30 new comets per year which ap¬ 
proach within 1 au of the Sun. This is about one order of magnitude 
larger than the discovery rate of 2.1 new comets (with q < 1 au) 
per year d Pones et alj|2004b[ ). If the PLPW model is correct and 
the planets do not change the flux of new LPC over a 10 Myr time 
span, and if most dynamically new comets with q < 1 au are in¬ 
deed discovered, then our simulations suggest a smaller population 
of OC comets, of order 10 11 comets, and thus a smaller total mass 
of the OC. 

We can also investigate the nonlinear effect caused by the 
combined perturbations from the tide and encounters (the so-called 
“syne rgy effect”) in delivering LPCs from the OC dRickman et alj 
l2008h . The synergy effect is measured by the fractional difference 
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Table 2. The number of LPCs injected into the LC (TV) and the number of dynamically new LPCs (n) for simulations of 10 million OC particles over the 
past/future lOMyr for the two encounter sequences. The subscripts c, t, and e refer to combined, tide-only, and encounter-only respectively. The last two 
columns give the two synergy parameters, s\ c and s new - 


Sequence 

Time span 

N c 

N t 

N e 

n c 

n t 

n e 

Ne/Nt 

rie/nt 

S1 C (%) 

Snew(%) 

Reliable 

-10 to OMyr 

19735 

18672 

1571 

7089 

6634 

281 

0.084 

0.042 

-2.3 

2.5 

0 to 10 Myr 

20151 

18873 

1560 

7147 

6708 

200 

0.083 

0.030 

-1.4 

3.3 

Basic 

-10 to OMyr 

19764 

18672 

1604 

7125 

6634 

279 

0.086 

0.042 

-2.6 

3.0 

0 to 10 Myr 

24571 

18873 

7080 

8441 

6708 

1525 

0.38 

0.23 

-5.6 

2.5 
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Figure 4. The relation between encounter-induced LPC flux and the two 
proxies ,7 (upper) and g (lower). For both panels, the red line denotes the 
optimized fitting function. The legends show the fitted functions of N e ( 7 ) 
and N e (g) and the Pearson’s chi-square for them (with 61 degrees of free¬ 
dom). 

in number of LPCs injected into the LC by the combined vs. en¬ 
counter and tide alone, i.e. 

sic = ( N c -Ne- Nt)/N c . (6) 

The synergy parameter for just the dynamically new comets is sim¬ 
ilarly 

Snew ~ (ttc ' Tit ')/ 7l c . ( 7 ) 

Values of s for both sequences and time spans are shown in Table[2] 
The synergy effect is not significant - only a few percent - for both 
sequences for both dynamically new and LC LPCs. This is small 
compared to the synergy effect found in the late stage of the 5 Gyr 
simulations for simulated encounters performed bv lRickman et alj 
|2008[) . This is probably a consequence of the relatively short time 
span of our simulations which does not permit us to include the 
long-term effect of massive encounters, plus the fact that we have 


very few massive encounters in our real data set iFouchard et alj 

ra . 

Generally, the synergy effect is negative for LC LPCs, but is 
positive for dynamically new LPCs. An encounter may perturb OC 
comets in the opposite direction to the tide, thus reducing the LC 
LPC flux compared to when either acts alone. Thus the negative 
synergy effect for the sequence may be caused by one or a few 
strong encounters. In contrast, the tide is not sufficient to inject 
comets into the observable zone directly (see the top right panel 
of Figure 0. Thus new LPCs tend to be injected by the combined 
impulse from the tide and encounters, leading to a positive synergy 
effect. 

We may also examine the flux ratio of encounter-induced 
LPCs over time and the evolution of the perihelia of some exam¬ 
ple LPCs. This is shown in Figure[5]for the reliable sequence. The 
shape of the flux ratio in the left two panels is driven by the obser¬ 
vational selection effects in our sample: there are fewer encounters 
the further we go into the past or future and so fewer injections 
at times far from the present. Nonetheless, injections can occur 
well after the encounter perihelion. According to the left two pan¬ 
els, the contribution of past encounters to the current LPC flux is 
about 6%, which is only a lower limit due to the incompleteness. 
The past-future asymmetry in the flux ratio and the sharp decrease 
around t = OMyr are caused by the time required by LPCs to 
reach their perihelia. Most encounters have a perihelion time within 
5 Myr from now, and have orbital period of around 5 Myr. Thus the 
flux ratio reaches its maximum at around 4 Myr and -1 Myr for the 
future and past sequences of encounters respectively. We see in the 
right panels that encounters can always intensify the variation of the 
perihelion of an LPC, and thus make it more likely to enter into the 
observable zone from outside of the LC within one revolution. This 
again explains why the synergy effect is positive for new LPCs. 

4.4 Angular distribution of the perihelia of LPCs and 
encounters 

Many previous studie s found that the dynamically new LPCs have 
anisotropic perihelia dDelsemmell 1 9871 : iMatese & Whitmirell201 ll : 
FBJ14). Here we reexamine this by investigating the effect of en¬ 
counters on the angular distribution of the LPC perihelia. 

The left panel of Figure [6] shows, in Galactic coordinates, the 
positions of the aphelia of encounter-induced LPCs (red crosses) 
and the perihelia of the stellar encounters (blue crosses) for the re¬ 
liable sequence of past encounters. We observe an East-West asym¬ 
metry in the directions of the aphelia of LPCs. This is due to the in¬ 
homogeneity in the distribution of strong encounters. We also see a 
significant difference between the locations of the aphelia of the ob¬ 
served LPCs (black points) and the simulated LPCs (red crosses). 
We further investigated the aphelia of LPCs induced by the tide and 
find a similar discrepancy. Given that the set of encounters we have 
is incomplete, and only accounts for about 5% of the LPC flux, it 
is not surprising that we see little correlation. The observed set of 
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Figure 6. An Aitoff projection in Galactic coordinates of the aphelia/perihelia of LPCs/encounters for the reliable sequence of past encounters (left panel) and 
HIP 103738 (Gamma Microscopii, right panel). In the left panel, the blue crosses show the perihelion of the encounters. The size of the cross is proportional 
to the number of LPCs the encounter has induced according to Table [T|( although a lower limit of the size is set for very weak encounters). The red crosses 
show the aphelia of dynamically new LPCs in the simulation and the black points show the aphelia of observed dynamically new LPCs. The black star symbol 
indicates the solar apex, and the dashed line shows the preferred tidal plane, the normal to which is the solar apex (FBJ14). The right panel is similar, but is 
only for HIP 103738, and does not show the observed LPCs and additionally shows the aphelia of LC LPCs in the simulation (yellow open circles). The blue 
curve is actually a sequence of points showing the encounter trajectory when its pericentric distance is less than 4 pc. Each point gives the star’s position at a 
equally-separated time interval of 0.001 Myr. 




t/Myr 


t/Myr 




Figure 5. The flux ratio, /, as a function of time (left panels) and the change 
of the perihelion of an example LPC (right panels) for the reliable sequence 
for future (top) and past (bottom) encounters. The left panels are similar 
to Figure [I] Note that f(t < 0) is set to be zero for future encounters. 
The right panels show the perihelion variation of an example LPC for both 
the tide-only (thin) and the combined perturbations (thick). The perihelion 
times of the five strongest future/past encounters are shown with vertical 
dotted lines and denoted by their Hipparcos IDs. 


LPCs is also likely to be incomplete (as what we see “now” only 
covers a few tens of years). 

We also see that the observed LPCs’ aphelia are located 
around the preferred tidal plane, where perturbations imposed on 
the OC are stronger due to the solar apex motion (FBJ14). In con¬ 
trast, the perihelia of encounters are only weakly associated with 
the tidal plane. A larger sample of encounters is required to con¬ 
firm the existence of the tidal plane in order to reduce the sample 


noise. We get similar results when we analyze the LPCs induced by 
the reliable sequence of future encounters. 

The right panel of Figure[ 6 ]shows the aphelia of LPCs induced 
just by HIP 103738 (Gamma Microscopii), one of the strongest re¬ 
cent past encounters. We find that the aphelia of both LC LPCs and 
dynamically new LPCs cluster near the trajectory of the encounter, 
although the geometry of the OC and characteristics of encounters 
complicates the spatial distribution of LPCs’ perihelia. In partic¬ 
ular, the aphelia of LPCs tend to concentrate around the location 
where the encounter was close to its perihelion. This is caused by 
the fact that an encounter generally has stronger geometric prefer¬ 
ence in imposing perturbations on the OC when it approaches its 
perihelion. 


5 DISCUSSION AND CONCLUSION 

We have conducted simulations of the perturbation of the Oort 
cloud in order to estimate the significance of known encounters 
in generating long period comets. According to simulations both of 
individual encounters and encounter sequences, the past encounters 
in our sample contribute about 5% of the current observed LPCs. 
Most LPCs were instead injected from the OC by the Galactic tide. 
We find that HIP 89825 (GL 710) and HIP 14473 are the strongest 
future and past encounters in our sample respectively, when we ex¬ 
clude objects with questionable data. We have also quantified the 
effect of other strong encounters. For most encounters, the peak of 
the encounter-induced LPC flux occurs about 2 Myr after the time 
of their perihelia. 

Consistent wit h the analytical estima tion s of the strength of 
encou nters given by Fouc hard et alJ 1201 lh and lMatese & Lissaueil 
d 2002 l) . our simulations show that the proxy g = M enc /(dp h u p h) 
can approximate the LPC flux induced by medium encounters rea¬ 
sonably well. Although some stars came (or will come) very close 
to the solar system, massive stars generally play a much larger role 
in perturbing the OC and increasing the flux of LPCs. In particu¬ 
lar, the recently noted close encounter by W0720 IMamaiek et alJ 
d2015h was extremely weak. 

We find that the synergy effect between encounters and the 
Galactic tide is not significant over 10 Myr time scales, although it 
could be significant over longer time scales. A comparison of the 
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simulated and observed flux of new LPCs suggests a population 
size of only 10 11 comets in the Oort cloud, which would require a 
smaller initial mass of the solar nebula than is often assumed. This 
may resolve a possible in consistency between the OC mass an d 
planet formation theories I Mor bidellil 120051 : [Duncan et alj l2008h . 
We also find that the perihelia of encounter-induced LPCs generally 
follow the trajectories of strong encounters. Since the tide cannot 
fully explain the anisotropic perihelia of LPCs, we expect that there 
are strong stellar encounters which have not yet been discovered. 
These may well be found by future astrometric surveys. 

Without accounting for the perturbations from the planets and 
observation bias, our model cannot precisely predict the current 
LPC flux. Yet our results on the flux of just dynamically new LPCs 
are more reliable because new L PCs have enough energy to pen¬ 
etrate the Jupiter-Saturn barrier IHeisler & Tremaine! 119861) . Al¬ 
though the new LPCs can al so form through the Kaib-Quinn mech¬ 
anism jKaib & OuinrJ[20Q9 f). the DLDW model has already taken 
this into account, and the planets can barely influence the flux of 
new LPCs over a lOMyr time span. 

The catalogs of encounters we have used are mainly derived 
from Hipparcos, limiting us to stars with V < 12 and with signifi¬ 
cant incompleteness. The Hipparcos catalogue has probably missed 
about one half of the encounters iFengiBoiS ). which have encoun¬ 
tered (or will encounter) the solar system over the past (or future) 
lOMyr. Gaia dPerrvman et al]|200ll) should detect more than 90% 
encounters over the past/future lOMyr and provide accurate as¬ 
trometric and spectrophotometric measurements for them dFenj 
120151) . Our work shows that future studies of encounters should 
focus not only on close encounters but also on massive and slow 
encounters. 
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